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Abstract 

We present a determination of the charm-quark mass in the MS scheme using the data combina- 
tion of charm production cross section measurements in deep -inelastic scattering at HERA. The 
framework of global analyses of the proton structure accounts for all correlations of the charm- 
quark mass with the other non-perturbative parameters, most importantly the gluon distribution 
function in the proton and the strong coupling constant as{Mz). We obtain at next-to-leading 
order in QCD the value mc{mc) = 1.15 ±0.04(exp) iQ QQ(scale) GeV and at approximate next-to- 
next-to-leading order mc{mc) = 1.24 ±0.03(exp) +j5;J]^ (scale) IJ!}? (theory) GeV with an accuracy 
competitive with other methods. 



'Permanent address: DESY, Notkestrafie 85, D-22607 Hamburg, Germany 



1 Introduction 



The charm quark is the lightest of the heavy quarks. Yet, the value of its mass nic is much larger 
than the scale Aqcd of Quantum Chromodynamics (QCD), i.e., iric ^ Aqcd- Thus, scattering pro- 
cesses involving charm quarks are subject to QCD dynamics at scales of the order of rric, where 
perturbative QCD predictions apply. This offers the opportunity to extract /n^. by comparing ex- 
perimental data for an appropriate observable to quark mass dependent theoretical predictions in 
perturbative QCD. This procedure does require some care, though. After all, quark masses are 
formal parameters of the QCD Lagrangian, but do not belong to the set of observables in Quantum 
Field Theory. Quarks and gluons do not belong to the asymptotic states at ? — )■ ±00 and already 
on grounds of the LSZ-theorem their mass is not the usual mass of a stable elementary particle, 
like the electron. No free quarks are observed in nature. As known from pertubation theory, the 
QCD corrections to the quark masses are renormalization scheme-dependent. Any quantification 
of these formal parameters necessarily assumes a definite choice of scheme a priori. 

In the past, high precision cross-section data from e+e^-coUisions have been the basis for such 
charm quark mass determinations. The available data from e+e -annihilation into hadrons span 
a large range of center-of-mass energies and can be used in QCD sum rule analyses based on 
perturbative QCD predictions to high orders in the coupling constant resulting in precise Mc values 
from scattering processes with time-like kinematics, see, e.g., Q. The recently available high 
precision data for charm quark production in deep-inelastic scattering (DIS) at the HERA collider 
now provide the attractive opportunity for a iric extraction from scattering processes with space-like 
kinematics. This is interesting per se for consistency tests of the Standard Model. Moreover, the 
precision now reached by the DIS measurements allows for an determination with an accuracy 
comparable to the one achieved in QCD sum rule analyses. 

In the present paper we use the new data combination of charm production cross section mea- 
surements in DIS at HERA [|2|| to determine m^. in the MS scheme by comparing to QCD predic- 
tions at next-to-leading (NLO) and next-to-next-to-leading (NNLO) order. We apply the formalism 
developed in Ref. fS\ and fit /n^ to the cross section data together with all other non-perturbative 
parameters, of which the gluon distribution function in the proton and the strong coupling constant 
OCi(Mz) in particular exhibit a significant correlation with m^-. For this purpose, we update the 
parton distribution function (PDF) analysis ABM 11 [4] with the new combined HERA data |l2l| 
included. Like ABMll, also the new variant of the present paper uses the MS renormalization 
scheme for (Xs{Mz) and the heavy-quark masses. It is performed in the so-called fixed-flavor num- 
ber (FEN) scheme for nf = 3 light quarks to be dealt with massless. The latter feature is rather 
important because in a global fit such as ABMll already the data for completely inclusive DIS 
measurements from HERA put significant constraints on the value of m^. due to the correlations 
mentioned. The FEN scheme allows for a well-defined description of open charm production in 
QCD, and the radiative corrections, i.e., the Wilson coefficients of the hard scattering process, are 
available exactly to NLO [l5]-|71 (see also Ref. [8j]) and to NNLO in an approximation for the most 
important parts, that is the gluon and quark pure-singlet Wilson coefficients lUl. 

The present study complements a previous determination of the c-quark mass [ITOl in the MS 
scheme based on data from the HI collaboration [[m[T2]| for open charm production. Those data 
are available in differential distributions so that the effect of value of on the extrapolation to 
the unmeasured region for the inclusive cross section has been carefully examined. In this way, 
Ref. IfTOl has obtained the MS mass md/Ur = mc) = mdmc) for the renormalization scale choice 
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Hr = mc at NLO to mdmc) = 1.27 ± 0.05(exp)+Q ^^(scale) GeV and at approximate NNLO to 
mdmc) = 1.36±0.04(exp)+[j;jJJ(scale) ±0.1 (theory) GeV, respectively. 

The present paper is organized as follows. In Sec. |2]we briefly recount the essential features of 
the data combination of Ref. |l2l. Sec.[3]contains the analysis and the new result for mc(mc) together 
with a detailed discussion on the impact of the new data set on the fit and the correlations of 
with the gluon distribution and the strong coupling as{Mz). We conclude in Sec. |4] emphasizing 
that the accuracy of the determination from DIS data becomes competitive with other methods, 
e.g., QCD sum rule analyses. 

2 Data 

The c-quark mass determination is conducted within the framework of a global analyses provided 
by the ABMl 1 fit [4J. The ABMl 1 analyses [4] has evolved from the previous ABKM09 fit 
and is based on world data for deep-inelastic scattering from HERA, and fixed target experiments 
and the Tevatron results on the Drell-Yan process. These data are supplemented by the recently 
published combined charm production cross sections in DIS at HERA Q and are used as input for 
the QCD analysis at NLO and NNLO. 

Reduced cross sections for charm production were measured in the kinematic range of pho- 
ton virtuality 2.5 < < 2000 GeV^ and Bjorken scaling variable 3 • 10"^ < Jc < 5 ■ 10"^ 
measurement was based on the combination of results obtained by using different charm tagging 
techniques: the reconstruction of D or D* mesons, the identification of muons from semi-leptonic 
decays of charmed hadrons or by exploiting the long lifetime in charmed hadron decays. The 
individual measurements were performed in different experimentally accessible (visible) phase 
space regions, depending on the experimental technique applied or on the different acceptances 
of the detector components used. For D-meson and muon production, the visible cross section 
measurements were extrapolated to the full phase space using predictions from perturbative QCD 
to NLO in the FEN scheme |l5l[6l[8]|. The quoted uncertainties in the extrapolation include those 
due to the variations of the factorization and renormalization scales, /jf, /Jr, simultaneously by a 
factor of 1/2 and 2 around the nominal scale, as well as of the charm quark mass in the range 
1.35 < m^^^ < 1.65 GeV for the pole mass m^^^ used in Refs. ^EM- 

The correlated systematic uncertainties and the normalization of the different measurements 
were accounted for in the combination procedure such that one consistent data set has been ob- 
tained. Since different experimental techniques of charm tagging were employed, the combination 
led to a significant reduction of the statistical and systematic uncertainties. However, due the com- 
bination procedure the information about the extrapolation factors and their uncertainties for the 
individual input data sets cannot be provided. Therefore, a detailed analysis similar to the pre- 
vious one performed in Ref. fTO\ taking into account the dependence of the extrapolation factor 
on the assumption of the charm mass in the underlying theory is not possible here. Instead, the 
results on the combined reduced charm cross sections at particular kinematical points (x, Q^) are 
used in the current analysis with account of the correlations of the uncertainties as provided by the 
experiments [[T4l . 
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3 Analysis 



The theoretical framework applied in the present analysis of the combined HERA data ^ essen- 
tially coincides with the one used earlier in the determination of the c-quark mass [[TOl with the 
HI data on open charm production [[ITl[T2l. We compute the heavy-quark contribution to the DIS 
cross section in the scheme with = 3 massless flavors in the initial state. The running-mass def- 
inition is employed for the heavy-quark Wilson coefficients, which comprise the NLO terms [|3l 
derived from the calculations performed with the pole mass definition [5 J and the NNLO terms [|9l. 
The latter are denoted by NNLOapprox in th^ following, obtained by interpolation between existing 
soft-gluon threshold resummation results and approximate relations for the Wilson coefficients at 
^ taking advantage of selected Mellin moments for the massive operator-matrix elements 
at NNLO given in Refs. lfT5] - [T9l . and the massless 3-loop Wilson coefficients EOll . The residual 
interpolation uncertainty which appears due to the finite number of Mellin moments being known 

is quantified by two options, c^^'^ and c^^'^ , for the constant terms in the Wilson coefficients 
at NNLO [9]. In the present analysis the shape of the NNLO correction is defined as a linear 
interpolation between these options using the ansatz 

cf) = (1 - Jyv)cf^'^ + ^^ivcf^'^. (1) 

The Fortran code OPENQCDRAD for the numerical computation of all cross sections in the present 
analysis is publicly available lf2T]| . 

Our determination of ruc is based on the 3-flavor ABMl 1 PDFs [!4|. However, those PDFs were 
obtained at the fixed value of mc{mc) = 1 .27 GeV. In order to provide a consistent treatment of the 
PDF dependence on iric we employ in the present analysis a set of /nc-dependent PDFs produced 
by interpolating between the variants of the ABMl 1 fit with the value of mc{mc) scanned over the 
range of 0.9 — 1.35 GeV. By fitting to the combined HERA charm data in this way we obtain the 
following c-quark mass values in the MS scheme 

mcK) = 1.20±0.05(exp) NLO, (2) 

mc{m,) = 1.30±0.04(exp) NNLOapprox- (3) 

Here the NNLO value corresponds to = —0.4 which provides the best agreement with the data 
in line with the approach of Ref. IfTOl . The experimental uncertainties in mc{mc-) are calculated 
by propagation of the errors in the data, taking into account the systematic error correlations. For 
the combined HERA data [2] they stem from 48 sources including the extrapolation of the visible 
charm production cross section to the full phase spaceH. This extrapolation is sensitive to the 
calculation details, such as fragmentation-model parameters, the PDFs, the value of m^, etc. The 
corresponding systematic errors encode the impact of the sensible variation of these parameters on 
the cross section values. Ideally, the extrapolation correction has to be calculated in the analysis 
iteratively, in parallel with fitting and the PDFs, as it has been done in the earlier determination 
of nic in Ref. [10] based on the selected set of the HI open charm production data. As discussed in 
Sec. [21 this approach is inapplicable in the present analysis because the necessary information about 
the visible phase space is lost in the combination of the HI and ZEUS data. The extracted value 

^The combined HERA data on open charm production with their systematic uncertainties used in the present 
analysis are available from http: //arxiv.org , as an attachment to the arXiv version of the present paper 
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of nic thus faces a procedural bias due to the fact that the extrapolation corrections are calculated 
for a fixed value of /n^. However, the corresponding uncertainty was estimated in Ref. [|2l by 
a conservative variation of the input used in the extrapolation correction. Therefore, the quoted 
experimental uncertainties in must exceed this bias. 

The central values of ntc in Eqs. (O and ([3]) are lower than those in our earlier determination in 
Ref. lUOl . In particular, this difference can be explained by a shift of the data obtained by the HI 
and ZEUS experiments in the process of their combination, cf. Ref. [2| for details. Besides, the 
NNLO correction employed in Ref. [|10| corresponds to the interpolation parameter = —0.6, 
which is somewhat different from the one obtained in the present analysis, which causes an addi- 
tional shift of mc{mc) in NNLO. However, in any case the values of mc{mc) in Eqs. ^ and ([3]) are 
compatible with the results of Ref. ifTOl within the uncertainties. 

To study the sensitivity of the determination to the particular choice of PDFs we repeat our 
analysis considering other 3-flavor PDFs. For this purpose, we take in all cases the nominal PDFs 
obtained with the fixed values of the c-quark mass. The NNLO values of iric obtained in this way 
demonstrate good agreement, cf. Tab. [T] At NLO only the ABM 11 [4] and GJR |l22l|23]| results 
coincide, while lower values are obtained in case of MSTW08 [24] and NN2 1 [|25l . This difference 
may partly appear due to a spread in the c-quark mass taken in different PDF fits. However, the 
difference between the ABM 11 results obtained with and without taking into account the Mc- 
dependence of PDFs is of 0(10) MeV, cf. Tab. [Hand Eqs. ^ and ([3]). This may point to other 
reasons for this difference. In fact, it is also correlated with the scheme used in the PDF fits. While 
the ABMl 1 and JR PDFs are based on the 3-fiavor scheme, the MSTW and NNPDF analysis are 
performed with different versions of a general-mass variable- flavor-number (GMVFN) scheme. In 
particular, this explains the difference at NLO between the MSTW and ABM/GJR results since 
the GMVFN scheme commonly deviates at NLO from the 3-flavor one to a larger extent than at 
NNLO. Recall also, that all PDF fits except ABMl 1, refer to the on-shell scheme for heavy quarks 
and compare to theoretical predictions using the pole mass ni^°^^. 





ABMii a 


JRfGJR) [1221231 


MSTW08 |24| 


NN21 [251 


NLO 


1.21 


1.21 


1.12 


1.01 


NNLO 


1.28 


1.27 


1.29 





Table 1: The value of mc{mc) in GeV obtained from the analysis of the combined HERA data on open 
charm production [[2] with different 3-flavor PDFs in NLO and NNLO. Note, the ABM 11 values are different 
from the ones in Eqs. ^ and © since the latter were obtained within the /Wf -dependent variant of the 
ABM 11 PDFs. 

Although Eqs. Q and ([3]) for mc{mc) are based on a consistent treatment of the PDF's c-quark 
mass dependence, the constraints on the variation of those PDFs with imposed by the data 
included in the ABMl 1 fit are not yet taken into account in the the determination of those numbers. 
To take advantage of the sensitivity of charm production in neutrino-nucleon scattering [[26 [[271 and 
the inclusive DIS to the charm mass we also perform the NLO and NNLO variants of ABMl 1 fit, 
which includes those data together with the HERA charm data of Ref. [^ added and the value of 
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mc{mc) considered as a fitted parameteiE]. From these versions of the fit we obtain the values of 

m,.(m,.) = 1.15 ±0.04(exp)+j!;jjJ(scale) NLO, (4) 

m,(m,) = 1.24±0.03(exp)l!!;!;3(scale)l!!:!!0(th), NNLOapprox, (5) 

where the NNLO value corresponds to dj^ = —0. 1 . This provides the best description of data with 
normalized by the number of data points (NDP), y}/NDP = 3459/3080 for the whole data set 
and y} /NDP = 61/52 for the combined HERA charm data, cf. also Fig.[T] At the same time, the 
option B of the massive NNLO correction of Ref. \9\ corresponding to Ja? = 1 is clearly disfavored 
by the data giving J^/NDP = 11 5/52 for the HERA charm data and J^/NDP = 3547 /3080 for the 
whole data set. Therefore we estimate the uncertainty due to the massive NNLO correction choice 
as variation between the values of mc(mc) obtained with dt^ = —0.1 and d/^ = 0.5 in Eq. ©• This 
yields the value of 0.07 GeV quoted in Eq. (|5]) as an estimate of the theoretical uncertainty. The 
scale uncertainty in mc(mc) is calculated as a variation due to a change in the factorization scale 
by a factor of 1/2 and 2 around the nominal value of ^/m]-. + KQ^, where K = 4 for neutral-current 
and K = 1 for charge-current heavy-quark production, respectively. For the NLO case both these 
variations lead to an increase in mc{mc) and we select the bigger shift as the uncertainty due to 
the scale variation. The NNLO scale uncertainty in mc(mc) is asymmetric and smaller than the 
NLO one, in line with the estimates of Ref. The experimental error of mc is reduced due to 
the constraints on the PDFs by the inclusive DIS data. The theoretical error due to missing higher 
order corrections is the dominant source of uncertainty in nic. 

The central value of mc{mc) obtained at NLO in Ref. |l2l for the combined HERA data including 
the data on charm production, i.e., mc{mc) = 1.26 ±0.05(exp) GeV, tums out to be bigger than 
our NLO result. It is important to note that this value is obtained from a scan of mc(mc) and not 
in a simultaneous fit of the PDFs and the charm quark mass. Also, partially the difference to our 
result can be explained by the different cuts on imposed in the analysis of Ref. |l2l and ours. 
By changing our cut of > 2.5 GeV^ to the cut of > 3.5 GeV^ used in (21 we get a shift of 
+0.03 GeV our NLO value of mc{mc) in Eq. dH). Another source for the difference is the data 
on dimuon production in neutrino-nucleon DIS [|26ll27l included in ABMll. By excluding this 
data set, we obtain a shift of +0.04 GeV for the c-quark mass in Eq. dH). Note, that the value of 
mc{mc) of Ref. [|3 is also systematically bigger than the NLO entries in Tab. [U which are obtained 
with fixed PDFs. Therefore the remaining difference between Eq. ^ and mc{mc) of Ref. ^ is 
evidently also related to particularities of the shape of HERA PDFs used in Ref. lO. 

Let us finally discuss a number of cross checks. Operating in the framework of a global analysis 
of the proton structure as provided by ABMll offers the possibility to account consistently for 
all correlations of the c-quark mass with non-perturbative parameters of the fit of which the gluon 
distribution function and the strong coupling constant a^(Mz) exhibit the strongest correlation with 
nic. We observe, that the shape of the gluon distribution obtained in the present fit is somewhat 
modified with respect to the ABMl 1 PDFs, cf. Fig. [2l However, the changes are basically found 
to be within the PDFs uncertainties. The sea distribution is affected to a lesser extent and the other 
PDFs are practically unchanged. 

^To allow for a variation of the factorization scale in the present analysis a cut of > 2 GeV^ is imposed on 
the data for dimuon production in neutrino-nucleon DIS Il26ll27l . while in the analysis of Ref. ||4| these data with 
~ 1 GeV^ were used. 
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The correlation of the fitted value of rric with the strong coupling constant as{Mz) is shown in 
Fig. [3]for a variation of the value of a^(Mz) in the range as{Mz) = 0. 1 10 - 0. 122. Recall, that the 
analysis of ABMll d has obtained a,(Mz) = 0.1180±0.0012 at NLO and a,(Mz) = 0.1134± 
0.0011 at NNLO as best fits. Fig. [3] demonstrates a remarkable stability of the c-quark mass both 
at NLO and NNLO. Considering a variation of 0.115 < av(Mz) < 0.119 the shift of Amdmc) is 
confined within an interval of 20 MeV for the NLO case and for a range of 0. 1 1 < (Mz) < 0. 114 
at NNLO within an interval of 10 MeV only. This is to be compared with the as{Mz) dependence 
inherent in QCD sum rule analyses. For example, for a variation of 0.113 < as{Mz) < 0.119 
Ref. [29] observes a linear growth of the value of mc{mc) with a maximal shift of l^dmc) — 
25 MeV (cf. Fig. 11a in [|29ll ). In contrast, the numbers for mc{mc) determined in Eqs. dU and ^ 
do not carry such bias with respect to the value of the strong coupling constant. 

To conclude the discussion we also convert the values of mc{mc) in Eqs. (H)) and ^ to the 
on-shell scheme. Using the well-known relations for the scheme transformation as encoded in [|30l 
and the values for a^iMz) of ABMl 1 at NLO and NNLO, we obtain 

mP°^^ = 1.35 ±0.05(exp)lj5;j;5(scale) NLO, (6) 

mP°i« = 1.59±0.04(exp)l!;:!!4(scale)l0:00(th), NNLOapprox- (7) 

As to be expected, the numerical values for mc°^^ axe larger than the values given in Eqs. © 
and ([5]) and those positive corrections grow in size, i.e., the shift of the central value amount to 
Amc(mc) = 200 MeV at NLO and Amc{mc) = 350 MeV at NNLO. The increasing spread between 
the numbers in Eqs. and (|7]) can illustrate the poor perturbative convergence of the pole mass 
scheme which is particularly pronounced at the low scales relevant for DIS charm production. 

4 Conclusions 

The new combined HERA data for charm production cross section measurements in DIS allows 
for a precise determination of the charm-quark mass in the MS scheme by comparing to QCD 
theory predictions in the FEN scheme at NLO and NNLO. Embedding the data analysis in a global 
fit takes advantage of a well-established theory framework and, simultaneously accounts for all 
correlations with other non-perturbative parameters, of which the gluon PDF in the proton and the 
strong coupling constant a^ ^Mz) are most important and have been studied in detail. 

The effect of the HERA DIS charm data on the extraction of (mt ) has been demonstrated 
in Eqs. ([2]), Yet, the full potential for a precision determination of nic^nic) unfolds in a global 
fit due the additional constraints imposed by the inclusive HERA data and those from neutrino- 
nucleon DIS. Thus, the best values for the c-quark mass are rric (me) = 1.15 ±0.04(exp) ljj;5JJ(scale) 
GeV at NLO and mdmc) = 1 .24 ±0.03(exp) +j!;55^( scale) IJjJ? (theory) GeV at approximate NNLO, 
cf. Eqs. (HI) and ([5]), although the accuracy of the latter determination still suffers from missing in- 
formation on the three-loop Wilson coefficients for neutral current DIS heavy quark production at 
small-x and small values of Q^. This implies an additional theoretical uncertainty on mc{mc) esti- 
mated to be in the range —70 < Aitic < MeV. The obtained values in Eqs. (01) and ([5]) are compati- 
ble with the previous analysis of Ref. ifTOl and with the world average mc(mc) = 1.275±0.025 GeV 
as summarized by the particle data group [[T]- The accuracy of the determination is competitive 
with other approaches, e.g., from scattering reactions in time-like kinematics. 
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Figure 1 : The combined HERA data on the reduced cross section for the open charm production [|2] 
versus x at different values of in comparison with the result of the present analysis at NLO 
(dashed line) and NNLO (solid line). A variant of the fit based on the option (A+B)/2 of the 
NNLO Wilson coefficients of Ref. [9j, cf. Eq. ([T]), is displayed for comparison (dotted line). 
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Figure 2: The relative change in the NNLO gluon (left) and non-strange sea (right) distributions 
obtained in the present analysis with respect to the ABMll PDFs (solid lines). The relative un- 
certainties in the PDFs are displayed for comparison (shaded area: ABMll, dotted lines: present 
analysis). 
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Figure 3: The values of mdmc) obtained in the NLO and NNLO variants of the ABMl 1 fit with 
the combined HERA charm data [2] included and the value of as{Mz) fixed. The position of the 
star displays the result with the value of as{Mz) fitted BUl. 
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